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ABSTRACT 

The abundance of clusters at the present epoch and weak gravitational lensing shear both constrain 
roughly the same combination of the power spectrum normalization eg and matter energy density £Im- 
The cluster constraint further depends on the normalization of the mass-temperature relation. Therefore, 
combining the weak lensing and cluster abundance data can be used to accurately calibrate the mass- 
temperature relation. We discuss this approach and illustrate it using data from recent surveys. 

Subject headings: cosmology: theory - large-scale structure of universe 



1. INTRODUCTION 

The number density of galaxy clusters as a function 
of their mass, the mass function, and its evolution can 
provide a powerful probe of models of large-scale struc- 
ture. Historically the most important constraint com- 
ing from the present day abundance of rich clusters has 
been the normalization of the linear theory power spec- 



trum of mass density perturbations (e.g. Evrard 1989, 


Frenk et al. 199C|, Bond & Myers 1991, |Hcnry & Arnaui 


1991 


, Lilje 1992|, Oukbir & Blanchard 1992|, Bahcall & Cer 


1993 


, [White, Efstathiou & Frenk 1993 


, Viana & Liddk 


1996 


. |Viana & Liddle 1999, Henry 2000 


). The normaliza- 



trast on scales 8 ft -1 Mpc, with the abundance constraint 
forcing models to a thin region in the Qm-&8 plane. 

Since the mass, suitably defined, of a cluster is not di- 
rectly observable, one typically measures the abundance 
of clusters as a function of some other parameter which is 
used as a proxy for mass. Several options exist, but much 
attention has been focused recently on the X-ray temper- 
ature. Cosmological N-body simulations and observations 
suggest that X-ray temperat ure and mass are strongly 



correlated with little scatter ( 


Evrard, Metzler & Navarro 


1996, Bryan & Norman 19981 


Eke, Navarro & Frenk 1998, 


Horner, Mushotsky & Scharf 1999, [Nevalainen, Marke- 


vitch & Forman 2000]). How well simulations agree with 



need to be resolved. On the simulation side there are the 
usual issues of numerical resolution and difficulties with 
including all of the relevant physics. On the observational 
side instrumental effects can be important (especially for 
the older generation of X-ray facilities) in addition to the 
worrying lack of a method for estimating "the mass". In 
this respect it is worth noting that there are numerous dif- 
fering definitions of which "M" an d "T" are to be related 
in the M-T relation ( [White 2001] )! 

With current samples the dominant uncertainty in the 
normalization i n fact comes from the normalization of the 
M-T relation (|Ekc, Cole fc Frenk 1996|, |Viana fc Liddk 
1996|, Ponahuc fc Voit 1999] , |Hcnry 2000[ pic7paoli, Scott 
fc White 20021, IScljak 20021), Or phrased another way, the 



cluster abundance is a sensitive probe of the normalization 
of the M-T relation. 

The abundance of clusters is, of course, not the only way 



to constrain the cosmological parameters. In this regard 
it is interesting to note that weak gravitational lensing 
provides a constraint on a very similar combination of f^f 
and us- Therefore, the two constraints can be combined to 
check for consistency of our cosmological model, to provide 
a normalization for the M-T relation, to probe systematics 
in either method and/or to measure other parameters not 
as yet included in the standard treatments. 

While the cluster constraint comes primarily from scales 
of about R = 10 /i -1 Mpc, current weak lensing surveys 
constrain somewhat smaller scales. These surveys probe 
scales between roughly 1 and 10 arcmin, which for source 
galaxies located at z ~ 1 in a ACDM cosmology corre- 
sponds to 0.7 h' 1 Mpc < R < 7h~ x Mpc. Therefore, weak 
lensing probes slightly smaller scales than clusters. As 
lensing surveys push to larger scales the overlap will be- 
come even better. 

In this paper we argue that a natural application of com- 
bining the cluster abundance and weak lensing constraints 
is to calibrate the M-T r elation for galaxy clusters (see 
also |Hu fc Kravtsov 2002J ). In Sec. | we define the M-T 
relation and derive how cluster abundance constraints de- 
pend on S]m and o~g. In Sec. Swe illustrate how combining 
the two constraints can fix the normalization of the M-T 
relation using two recently obtained data sets. Finally, in 
Sec. |] we discuss this approach further. 

2. THE MASS-TEMPERATURE RELATION 

Throughout we shall be interested in the abundance of 
massive clusters at low redshifts, so we parameterize the 
M-T relation as 



M(T, z) 
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-3/2 



(1) 



where M15 = 10 15 M Q , A c is the mean overdensity in- 
side the virial radius in units of the critical density, which 
we compute using the spherical top-hat collapse model, 
and E 2 = f2 M (l + z) 3 + £1 A + £l k (l + z) 2 . T* is the normal- 
ization coefficient that we seek to constrain; it roughly cor- 
responds to the temperature of a M — 7.5 x 10 13 M Q 
cluster. If measured in keV, the value of T* is precisely 
equivalent to j3 from Pierpaoli, Scott & White (2002) and 
is 1.34/ T of Bryan & Norman (1998). 

Let us explore the sensitivity of cluster abundance on 



2 



£Im and o% . The Press-Schechter formula gives the num- 



ber o f collapsed ob jects dn per mass interval rflnM (Prcsi 
& Scr jechter 197^ ); we dehne M{M,z) = dn/dhxM. Fur- 
ther defining i/ = S c /a(M, z), where <j(M,z) is the rms 
density fluctuation on mass-scale M evaluated at redshift 
z using linear theory and S c rs 1.686 is the linear threshold 
overdensity for collapse, we have 



Af(M, z) 



2 pm d\na(M, z) 



7T M 



din / 



1/ exp {-v 2 /2) 



(2) 



where is the present-day matter density. Assuming 
we are dealing with t he current cluster abundance, z ~ 0. 
Following Pen (1998), for the mass scales of interest we 
can approximate a(M) oc M _a where a ~ 0.27 for the 
currently popular ACDM cosmology. 

Let us examine the dependence of Af on Qm , as and M. 
Ignoring the term d\ua / d\nv (which slowly varies) one 
obtains 



SAf sn M 



5 a s 



SM 
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<7g M 

(3) 

Setting the left-hand side to zero and using the fact that 
SM/M = —3/2 <5T*/T*, for our fiducial cosmology and 
massive clusters (M ~ 10 15 Ii~ 1 Mq, or v ~ 2) we haveR 



T* oc (o-gfiS/)- 1 - 1 . 



(4) 



Therefore, measurements of the cluster abundance at 
the present epoch constrain a degenerate combination of 
T* and og f^/- One °f them cannot be determined with- 
out knowing the other. Thankfully, weak lensing happens 
to measure roughly this combination of Qjm and as ac- 
curately, an d the orthogonal combin ation much less accu- 



rately (e.g. Bernardeau et al. 1997). Consequently, weak 



lensing in conjunction with cluster abundance can be used 
to constrain T» quite strongly. 

3. WEAK LENSING PLUS CLUSTERS: AN EXAMPLE 

As a more concrete example of these ideas, let us exam- 
ine what value of T* is required to bring current cluster 
and weak lensing results into agreement. This analysis will 
necessarily be illustrative, but is already quite enlighten- 
ing. 

3.1. The cluster data 

We compute as using a Monte-Carlo method follow- 
ing the steps outlined in Pierpaoli, Scott & White (2002). 
Since some of the details have changed we sketch the pro- 
cedure here. 

We use the HiFluGCS cluster sample of Reiprich & 
Bohringcr (1999), restricted to clusters with 0.03 < z < 
0.10. For simplicity we do not include 'additional' clus- 
ters of lower flux/temperature which could scatter into 
the sample. The cosmic microwave ba ckgro und (CMB) 
frame redshifts from Struble & Rood (1999) were used 
when available and so were the two-component tempera- 
tures published in Ikebe et al. (2002). For each Q m we sam- 
ple from a distribution of cosmological parameters includ- 
ing h, n and (the normalization of the M-T relation). 

1 Note that the dependence of M (or T«) on is stronger for 
more massive clusters; a mnrp detailed analysis gives T» oc cr -5 / 3 for 
the most massive clusters (Evrard et al. 2002). 



For each such realization we generate 50 mass functions, 
where the temperature is chosen from a Gaussian with the 
mean and variance appropriate to the observational value 
and errors, and a scatter of 15% in mass at fixed T is as- 
sumed for the M-T relation. Using the mean values of the 
M-T relation and the L-T relation from Ikebe et al. (2002) 



L x = 1.38 x 10 J 



/ kT \ 2 5 
V IkeV J 



h~ 2 W 



(5) 



we compute the volume to which clusters of mass M could 
be seen above the flux limit /n m = 1.99x 10~ 14 erg s _1 cm~ 2 
of the survey. For each realization of the mass function we 
compute the best fitting as by maximizing the Poisson 
likelihood of obtaining that set of masses from the theory 
with all parameters except as fixed. The mass function 
can be compute d usin g either the Press-S chech ter (1974), 
Sheth-Tormen ( 1999 ) or Jenkins et al. ( 2001 ) formulae. 
We have used the Sheth-Tormen prescription throughout, 
with the mass variance a 2 (M) comput ed usi ng the trans- 
fer function fits of Eisenstein & Hu (|1999| ) and masses 
converted from Afisnn. to Ma assuming an NFW profile 
( |Navarro, Frenk fc White 1997| ) with c = 5. The best fit- 
ting as is corrected from z to z — 0. The mean of the 
50, z = normalizations is then taken as the fit for that 
set of cosmological parameters (since the error from Pois- 
son sampling is completely sub-dominant to the error in 
the M-T normalization we do not keep track of it here). 
When quoting a best fit for a given triplet of (fl m , as, 
T»), we marginalize (average) over the other cosmological 
parameters h and n. 

3.2. The weak lensing data 

As an example of weak lensing measurements, we use 
shear measurements obtained using Keck and William 



Hcrschel telescopes (Bacon et al. 2002). These joint mea- 
surements used two independent telescopes covering 0.6 
and 1 square degrees respectively, and enabled careful as- 
sessment of instrument-specific systematics. The authors 
compute the shear correlation function, and compare with 
the theoretical prediction. Assuming the shape parameter 
T = 0.21, the results are well fit by 



08 



0.97 ±0.13 



(6) 



which captures the total 68% CL error: statistical, red- 
shift uncertainty and uncertainty in the ellipticity-shear 
conversion factor. These results are consistent with other 



recent measurements of cosmic shear ( van Wacrbcke et al 



20021. [Rcfrcgicr, Rhodes 
Gladdors 2()02|) ~ 



Groth 2002, Hoekstra, Yee fc 



3.3. Calibrating the M-T relation 

Fig. [l] shows the constraints in the £Im-&8 plane. The 
cluster constraint has been marginalized over h and n as 
explained above, and plotted for three different values of 
T*. We have checked that the allowed ranges for h and 
n are wide enough so that essentially all of the likelihood 
is contained within those ranges. The weak lensing con- 
straints assume the shape parameter T = 0.21. Note that 
the constraint regions from the two methods are indeed 
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Fig. 1. — 68% CL un certainty contains in the Q M -ag plane, for 
a weak lensing survey ( [Barnri et id ?Q0?| ) and 95% ,CL uncertain- 
ties for a cluster survey ( Keiprich & Bohringer 1999). The cluster 
results are shown for three different values of the mass-temperature 
normalization parameter T*, and marginalized over n and h. The 
degeneracy regions for the two methods are very similar, which in 
principle enables an accurate determination of T, . 



parallel, with very similar degeneracy directions. This en- 
ables an accurate determination of the normalization T* . 

In the example above, we see that a relatively low T* 
is preferred (T* < 1.7 keV) in order for cluster results to 
agree with the weak lensing results. While systematics 
in both methods could still be important, it is interesting 
to note that this result is in line with most earlier esti- 



mates ( Evrard. Mctzlcr fc Navarro 1996 ,_ Eke, Navarro & 



1998 , Bryan fc Norman 1998 , Yoshikawa, Jing & 



iOOO), while it disa grees with values adopted more 



Frenk 
Suto _ 

recently (e.g. |Seljak 2002[ ). 

The fact that cluster abundance and weak lensing probe 
different scales opens a possibility that one might be able 
to secure the agreement between the two methods by vary- 
ing the shape of the power spectrum or the spectral in- 
dex n rather than the M-T normalization. Unfortunately 
the constraints we have combined above have individually 
been marginalized over h and n. Ideally, one would com- 
bine the cluster and weak lensing likelihood functions and 
then marginalize over the relevant parameters to get the 
probability distribution of : 

P(1V) = J £ c i UB (T*,fiM,0- 8 ,n,/i) 

x Avl(^m, erg, n, h) (KIm das dndh. (7) 

Then the results would be manifestly independent of the 
power spectrum parameters. We do not have the ability 
to perform such an analysis here. 

Note, however, that the scales probed by lensing and 
clusters are quite close, separated an order of magnitude 
at most. For example, it would require a spectral tilt of 
n ~ 1.2 to make the recently obtained "low" normalization 
from cluster abundance (erg ~ 0.6) agree with the "high" 
normalization from weak lensing (erg ~ 0.9), and such a 
high value of n is already disfavored by recent CMB ex- 



periments flBalbi et al. 2000|, INetterfield et al. 2002j , [Pryke 
et al. 2002| , |Sicvers ct al. 2002| ). 



4. CONCLUSIONS 

There has been a lot of discussion recently regarding 
the value of cluster normalization er 8 . W h ile the "old" re- 
sults favor erg ~ 1 (Viana & Liddle 1999, Pierpaoli, Scott 
& White 2002 and references therein), several new clus- 



malization ( Rci 


jrich & Bohringer 199G, Borgani ct al 


2001, Viana ct al 


. 2002, Scljak 2002, |Ikebe et al. 2005 , Bah- 


call et al. 2002|). 

by the combiner 
and CMB data 


Thc lower normalization is also favorec 
analysis of 2dF Galaxy Redshift Survey 
'Lahav et al. 2001). On the other hand, 



recent weak lensing results (van Wacrbckc ct al. 2002 



Yee 



Ba 



con et al. 20021 |Refregier, Rhodes fc Groth 2002j , |Hoekstra, 



Gladdcrs 2002) tend to favor a higher value of a%. 
The cause is of this discrepancy between various measure- 
ments has not been identified yet; one candidate is larger 
than anticipated systematic errors in one or both methods. 
Another possibility is the bias in the relation between the 
mass and the observable quantity — temperature or lumi- 
nosity — used to construct the abundance of clusters. 
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Fig. 2. — Estimates of the M— T normalization T* collected from 
the literature. The red (light) points are estimates from simulations, 
while the black (dark) points are from the observations. Points with 
no error bars had none quoted. The shaded region is roughly our 
favored range of values of T» . 

The cluster abundance constraint on a$ crucially de- 
pends on the M-T normalization T*. Figure [2] summa- 
rizes the current status of our knowledge of T„ . It shows 
seven determinations from N-body simulations and three 
from direct observations, as compiled in Pierpaoli, Scott 
& White (2002) and Muanwong et al. (2002). The shaded 
region is roughly our favored range of values of T* . Points 
without error bars had none quoted, and the three ob- 
served values of T* assumed the isothermal-/? model. The 
measurement due to Muanwong et al. corresponds to their 
"radiative" and "preheating" cases that are cooling-flow 
corrected, while the value due to Pierpaoli, Scott & White 
is an average over the simulations. The large discrepancy 



4 



between the different measurements is apparent, and it 
also appears that the observed values are systematically 
higher than the ones obtained from simulations (see Muan- 
wong et al. 2002 for further discussion). 

We argue here that the cluster abundance - weak lens- 
ing complementarity can be used to cross-check the M— T 
relation. By combining recent weak lensing constraints 
from Bacon et al. and the HiFluGCS cluster sample of 
Reiprich & Bohringer, we have demonstrated the util- 
ity of this method. While potential systematic errors in 
both data sets are still a concern, the example we used 
prefers relatively low values of the M-T normalization 
(T* < 1.7 keV). We conclude that future weak lensing 
surveys (Vista, LSST, SNAP) combined with new cluster 
data from Chandra and XMM-Newton observations will 
provide a strong probe of the M-T relation. 
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